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Abstract 



An enhanced geothermal system (EGS) expands the potential of geothermal energy by enabling the exploita- 
q tion of regions that lack conventional hydrothermal resources. The EGS subsurface system is created by 
& engineering enhanced flow paths between injection and production wells. Hydraulic stimulation of existing 
W fracture networks has been successfully achieved for unconventional geothermal resources. More recently 
£N proposed concepts increase the use of drilled wellbores in hard rock to connect the injection and production 
wells. The present work investigates the long-term thermal effects of deviated geothermal wellbores and stud- 
ies how the cooling of the borehole wall results in thermally induced tensile fractures. The results show that 
^induced fractures are created by a combination of in situ and thermal stresses, and that the extent to which 
q thermally induced tensile wall fractures are created largely depends on how the wellbores are oriented with 
^respect to the pre-existing stresses of the reservoir. If the system is not optimized with respect to in situ 
^ stresses, the risk of wellbore instability becomes severe within less than a year of production. In contrast, 
O if the orientation of the wellbores is optimized, thermally induced instabilities can be completely excluded 
# ui as potential risks for the operational lifetime of the system. Furthermore, our results show that the thermal 
^failure process strongly depends on the temperature of the injected water but is only weakly affected by the 
Qninjection rate. 
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1. Introduction 

Enhanced geothermal systems (EGSs) have 
emerged as a promising new type of geothermal 
power technology. An EGS uses hydraulic fracturing 
to create a fracture network in low-permeability rock 
layers to enhance the connectivity between typically 
deviated injection and production wells ( |Abe et 
all [T999|). Still in the development stage, EGS 



of deviated wellbores to ensure connection between 
the production and injection wells (Fig. [TJ. One 
approach avoids the use of hydraulic stimulation, 
thereby reducing the risk of induced seismicity, by 
interconnecting the production and injection wells 
via drilled wellbores ( |Holmberg et al. , 2012). The 
other approach uses drilled boreholes to improve the 
connection between the wellbores and the existing 



technology has been criticized because of its po- 
tential to cause the initiation of larger seismic 
events (Maje r et aL} |2007| ). Recently, two alternative 



EGS concepts have been presented ( [Holmberg et 
aL} |2012[ |Zhang et~aL} |2012| ) that increase the use 
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or stimulated fracture networks (Zh ang et al4|2012| ). 
The present study refers to these new concepts as 
wellbore-based EGSs (WEGSs). The operational 
lifetimes of both EGSs and WEGSs largely depend 
on long-term effects, such as mineral scaling, cor- 
rosion, and fracturing resulting from thermal stress, 
which could alter system permeability. The present 
article focuses on how long-term thermal effects can 
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result in borehole instabilities due to fracturing. 




Figure 1 : Illustration of a WEGS concept based on connected 
wellbores. The injection and production wells are intercon- 
nected by drilled wellbores (heat exchangers) that are tilted by 
an angle of 6 with respect to the z axis. The length of each 
heat exchanger is 2000 m. Zb represents the z-direction in the 
borehole coordinate system, while z is the vertical direction in 
the geographic coordinate system. The maximum depth of the 
wells is 5000 m below the surface. The temperature at 5000 
m is assumed to be 160 °C, while the water enters the heat ex- 
changers with a temperature of T^j . 

Stress is induced in a reservoir that is exposed 
to compression or tension. The sources of inter- 
nal stress can be forces that are created due to a 
change of the reservoirs temperature or, in the case 
of a porous medium, a change in fluid pressure. 
Large pre-existing (in situ) stresses are also typi- 
cally present in the reservoir due to tectonic forces 



and the overburden pressure ( |Zoback[ |2007). Tensile 
or compressive rock failure is observed if the stress 
exceeds the rock strength (ZobackJ |2007| ). Model- 
ing this fracturing is complicated, because the pro- 
cess depends on several coupled mechanisms that act 



on different length and time scales ( |Hansen[ [2005 
Hals and Berre[ |2012| ). Fracturing has been stud- 
ied extensively over the last few decades with re- 



and Gonzalez} [T984t |Peska and Zoback} [1995) ). In 
the petroleum industry, one of the major operational 
concerns related to drilled wellbores has been the 



prevention of borehole wall failure (Cooper, 1994). 



This type of failure is initiated by pressure and tem- 
perature gradients that are created around wellbores 
during drilling (Per kins and Gonzalez} |1984| ). Sev- 



eral previously published studies have investigated 
the stability of inclined boreholes. In particular, 
pressure-induced failures have been the focus of oil- 
well studies. In addition to numerical and experi- 



mental studies (Daneshy, 


1973i Yew and Li, 1 1988; 


|Baumgartner et al., 1989^] 


Rawlings et al. , 


1993 


),the 



problem of borehole wall failure has been studied an 
alytically via linear elastic models ( Bradley} |1979 



|Aadnoy[[T987l ). |Peska and Zobackj ( [T995] ) presented 
a thorough investigation of borehole stabilities based 
on arbitrary borehole orientations in a wide variety of 
in situ stress states. They employed analytic expres- 
sions that describe the principal effective stresses at 
the wellbore wall of a cylindrical well while consid- 
ering fluid pressure variations ( |Hiramatsu and Oka 



1968 ; |Fairhurst[ [1968] ). Similar analytic expressions 
have been used in studies of thermally induced fail- 
ures ( [Myklestadl 1 1 942} [Perkins and Gonzalez[ 1 1 984 



Zoback[ |2007| ). The effects of thermal fracturing 



have also been of concern in carbon capture and 
storage (CCS) technology. Recent studies have in- 
dicated that thermally induced stresses significantly 
influence the safety of the seal and the fracturing 
of the caprock ( |Luo and Bryant} |2010[ |Goodarzi 
eTaL} |20T(j} |20T2| ). In the context of geothermal 



reservoirs, several studies have focused on thermally 
induced stresses, mainly concentrating on heat ex- 



traction enhancements due to thermal cracks (Mur 



phy} [T979} IDemuth and Harlowj [T980| ). Thus far, 



gard to the stabilities of deviated boreholes (Perkins 



few works have attempted thorough stability analy- 
ses of geothermal wells that investigate the interplay 
between in situ stresses, borehole orientations, and 
long-term thermal effects. 

The present paper combines numerical and ana- 
lytical methods in the first theoretical study of the 
long-term thermal effects of deviated geothermal 
wellbores in regions where convective fluxes be- 
tween geothermal reservoirs and wellbores are neg- 
ligible. Temperature evolution in the geothermal 
well and rock reservoir are simulated using a numer- 



ical COMSOL Multiphysics model ( |Littmarck and 
Saeidil[T986), while stresses at the wellbore wall are 



calculated from analytic expressions valid for cylin- 
drical wells. We show that the degree of thermal 
tensile fracturing is extremely sensitive to how the 
drilled wellbores are oriented with respect to in situ 
stress states. In the case of a large horizontal stress, 
the wellbores should be oriented along the maximum 
horizontal stress direction. This orientation nearly 
eliminates the risks for thermally induced borehole 
instabilities, even after several years of injecting cold 
water. In contrast, when the wellbores are oriented 
along the minimum stress direction, wellbore insta- 
bilities are observed within less than one year. In 
addition, we show that the degree of thermal fractur- 
ing strongly depends on the temperature of the in- 
jected water but is only weakly affected by the injec- 
tion rate. The strong fluid-temperature dependency 
of the thermal failure process allows one to control 
the degree of fracturing for a given in situ stress state. 
Our results provide important insights for optimizing 
WEGSs as well as oil and geothermal wells, in gen- 
eral. We expect that the results can be applied to CCS 
technology to enhance the safety of geological stor- 
age sites. 

This paper is organized as follows. In Section |2j 
we present the theory and governing equations that 
model the geothermal system. Section [3] provides a 
description of the model system under consideration, 
while Section]?] presents our findings. We conclude 
and summarize our results in Section [3 

2. Governing Equations 

The following sections, present the equations that 
govern the temperature evolution of the reservoir and 
the heat transfer process between the hot rock and 
cold water that circulates through the wellbores. In 
addition, the analytic expressions for the stresses that 
act on the wall of deviated, cylindrical wellbores in 
the presence of in situ stresses, fluid pressure, and 
thermal stresses are presented. 

2.7. Temperature Equations 

The temperature equations for the system consid- 
ered herein are derived from the conservation of en- 
ergy, which is described by two coupled processes: 



the energy change of the fluid and the energy change 
of the rock. These two processes are coupled by 
the heat transfer that occurs at the borehole wall. 
This heat transfer is a complicated process that de- 
pends on several different parameters, including the 
velocity and the temperature of the fluid, the min- 
eral composition of the rock, and the wellbore prop- 
erties. The system is characterized by the following 
two temperature equations that describe the tempera- 
ture, Tf(r, t), of the fluid inside the wellbore, Q w , and 
the temperature, 7V(r, t), of the rock in the region, Q r , 



outside the wellbore (Bejan and Kraus, 2003): 



PfCpJ [ih + V " V7> ) = V " (* /V7> ) ' X G Qw ' (1) 



p r C 



p,r 



~dt 



= V • (K r VT r ) , x G Q r . 



(2) 



Here, p f (p r ), c p j (c Ar ), and *y (K r ) are the density, 
isobaric specific heat capacity, and thermal conduc- 
tivity of the fluid (rock), respectively, while v is the 
fluid velocity. These two equations are coupled by 
the heat transfer between the fluid and rock. To low- 
est order in the temperature difference (T f - TV), the 
heat flux across the borehole wall is given by (|Bejan| 



and Kraus ( 2003): 



q = hiTooj - 7V)n, 



(3) 



where n is the surface normal of the wall pointing 
into the rock reservoir, h is the heat transfer coeffi- 
cient, TV is the temperature of the borehole wall, and 
j is the bulk temperature of the fluid in the well- 
bore. The heat transfer coefficient can be expressed 
in terms of the Nusselt number, Nu, of the fluid by 



the following relation (Bej an and Kraus[ [2003): 



Nu 



hD h 



(4) 



Here, D h is the hydraulic diameter, which is equal to 
the borehole diameter, D, for a cylindrical wellbore. 
Several empirical models for the Nusselt number ex- 
ist. In the present study, we calculate Nu from the 
Dittus-Boelter equation (see Section [3] for further de- 
tails). 
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2.2. Effective Stresses and Failure Criterion 

In real rock reservoirs, in situ stresses are present. 
The instability of a borehole largely depends on its 
orientation with respect to in situ principal stress di- 
rections. In the following analysis, we present an- 
alytic expressions for the effective stresses on the 
borehole wall for arbitrary borehole orientations and 
principal in situ stress directions. Further details can 



be found in Peska and Zoback (1995) and Zoback 



(2007), and the references therein. 



Let S ij denote the in situ stress tensor represented 
in the borehole coordinate system, where the z b axis 
points along the borehole axis (Fig.[T]), the x b axis is 
in the plane perpendicular to the borehole axis (ori- 
ented towards to the bottom side of the borehole), 
and the y b axis is in the same plane, but is perpendic- 
ular to x b (see |Peska and Zoback|fll995] ) for an illus- 
tration). The effective stress tensor in the borehole 
coordinate system is then, cr y = 5 y - bP p 5ij, where 
P p is the pore pressure, b is the Biot- Willis parame- 
ter, and 6ij is the Kronecker delta. A positive stress 
implies that the system is under compression. The 
effective stresses at the borehole wall of a cylindrical 
wellbore are: 

o~ ZbZb = o- 33 - 2v (cr n - 0-22) cos(2#) - 
4v(r 12 sin(20), 



O-ffi = (Tix + 0~22 ~ 2 ((Tn 

4cr 12 sin(2#) - Ap - 



cr 22 )cos(2§)- 



T~ 6zb = 2 (cr 23 cos(§) - o~ u sin(0)) , 



(Trr = A/7 + of . 



(5) 



Here, 6 is the angle (around the borehole) measured 
counterclockwise from the x b axis, r is the radial di- 
rection away from the center of the cylindrical well, 
Ap is the difference between the borehole fluid pres- 
sure and pore pressure in the rock, and cr^J and <x~~ r 
are the thermally induced stresses. In the present 
study, we are only concerned with the value of cr~~\ 
which can be reduced to the following expression in 



the steady state (Zoback, 2007) 



Here, v is the Poisson's ratio, a t is the linear coeffi- 
cient of thermal expansion, E is the Young's modu- 
lus, and AT(t) = T bwJ - T bw (t) is the difference be- 
tween the initial temperature, T bwJ , of the borehole 
wall and the time-dependent temperature, T bw (t), of 
the heated or cooled borehole wall. Although the 
above equations assume that Ap and AT are constant 
along the borehole axis, the expressions constitute 
good approximations for more general situations if 
the spatial variations of AT and Ap along the bore- 
hole axis are weak compared to the spatial variations 
along the radial direction. 

From Eq. @, we see that increasing the borehole 
fluid pressure or cooling the rock reduces the tangen- 
tial stress cr~Q~Q. For sufficiently large values of Ap and 
AT, the system may evolve from a compressive state 
into a tensile state. Because the tensile strength of 
rock is relatively weak, a common tensile failure cri- 
terion (for a fixed 6 direction) is when the minimal 
effective stress in the plane tangential to the borehole 
becomes negative ( |Zoback[ |2007| ). The minimal ef- 
fective stress is given by: 



(Tt 



°~ z b Zb + °~00 



^) 2 + 4r 2 



Ozbi 



0% = a t EAT/(l-v). 



(6) 



(7) 

Therefore, the present analysis assumes the creation 
of tensile wall fractures when cr t min < 0. 

3. Model Description 

This section provides a description of the model 
system used herein to study of thermally induced in- 
stabilities in geothermal wells. 

3. 1 . COM SOL Multiphysics Model 

A COMSOL Multiphysics (version 4.3) model is 
used to simulate the temperature evolution of a devi- 
ated wellbore cooled by injected fluid. In the case of 
multiple wellbores, we assume that the separations 
between the boreholes are sufficiently large to avoid 
thermal interaction, such that each borehole can be 
modeled separately. We only simulate the temper- 
ature evolution along a 2000m-long wellbore seg- 
ment. In the WEGS concept illustrated in Fig. [TJ 
such a wellbore segment is provided by one of the in- 
clined boreholes that connects the injection and pro- 
duction wells. The COMSOL model of the 2000m- 
long borehole segment is illustrated in Fig [2j The 



system is modeled in the borehole coordinate system, 
in which the location of the well segment is always 
between 3000 m < Zb < 5000 m, irrespective of the 
tilting angle, 6. The system is assumed to be rota- 
tional symmetric around the Zb axis, and the radius is 
50 m. Rotational symmetry is also a good approx- 
imation in the case of a tilted wellbore, 6 > 0, be- 
cause the difference between the temperature at the 
center of the system and at the outer surface (for a 
fixed Zb value) is ~ 50 sin(0)dT/dz. For the temper- 
ature gradient, dT /dz, considered here, this situation 
yields a negligible temperature difference (less than 



z =3000 



1.5 K). The diameter of the wellbore is 0.1 m ( [Holm- 
berg et al.[ |20T2t |Zhang et al.j [20T2| ). The outlet of 
the wellbore segment is always located z = 5000 m 
below the surface (in the geographic coordinate sys- 
tem), while the depth of the inlet depends on the ori- 
entation of the borehole and is related to the tilting 
angle by: z = 3000 + 2000(1 - cosfl) m. The tem- 
perature evolution of the fluid and rock matrix are 
described by Eq. ([I]) and Eq. ([2]), respectively. The 
injected fluid is water. The fluid properties, including 
the density, thermal conductivity, and heat capacity, 
are given by the built-in COMSOL functions for wa- 
ter. The fluid velocity, v, is assumed to be constant 
throughout the entire wellbore and is calculated from 
the injection rate and fluid density at a temperature of 
60 °C. The rock matrix has the following properties, 



which are typical for Fennoscandia (Slagstad et al 



2009): the heat capacity is c PJ = 850 J/kgK, the den- 
sity is p r - 2600 kg/m 3 , and the thermal conductivity 
is K r = 2.8 W/mK. The initial temperature, Ti(zb) 9 of 
the rock reservoir in Fig. [2] depends on how the well- 
bore is tilted with respect to the vertical z axis and is 
given by: 



(8) 



dT 

Ti(z b ) = T + -i—Zb, 
dzb 

dT 

T = T s + 5000(1 -cos(#))— , 

dz 

dT dT 

= — cos(0). 

dzb dz 



Here, 6 is the tilting angle, T s = 4° is the surface 
temperature, and dT/dz = 0.0312 K/m is the tem- 
perature gradient. The water in the wellbore is ini- 
tially assumed to be in thermal equilibrium with the 
rock (before the simulation starts and cold water is 



z =5000 




7(3000) 



I I I 



n-V( Kf T)=0 



T=160° 



50 m 



0.05 m 



Figure 2: Illustration of the COMSOL model. The system is 
rotational symmetric around the Zb axis. The model consists 
of a rock reservoir with a radius of 50 m. The drilled well- 
bore, which has a diameter of 0.1 m, is located at the center 
of the reservoir. The borehole and rock reservoir are sepa- 
rated by a thin, thermally resistive layer (red line). The sys- 
tem is thermally insulated at the upper and lower surfaces (i.e., 
n-V(/cyT) = 0, where n is the surface normal). The outer surface 
temperature is assumed to be the initial temperature, Ti(zb), of 
the reservoir. The injected water has a temperature of T^j and 
maintains a uniform velocity profile through the wellbore. 



injected into the wellbore with a temperature of r in j). 
The boundary conditions are shown in Fig. [2j The 
upper and lower surfaces of the rock matrix are ther- 
mally insulated; thus, it is assumed that the thermal 
currents along these two surfaces are oriented along 
the radial direction. To model a rock system that is 
in contact with an infinite rock reservoir, the temper- 
ature along the outer surface is set equal to the initial 
temperature of the reservoir. 

Heat transfer between the rock and fluid is mod- 
eled by a thin, thermally resistive layer at the bore- 
hole wall, such that the heat flux across the layer is 
given by: 

Tf-Tr 

q = K f 



d 



-n. 



(9) 



Tf and T r are the temperatures of the fluid and rock 
on each side of the thermally resistive layer, respec- 
tively, and d is the thickness of the layer. The thick- 
ness is determined by equating Eq. ([9]) and Eq. ([3]), 
which yields d = D/Nu. The Dittus-Boelter equa- 
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tion, known to be valid for fully developed, turbulent 
flow, is used to calculate the Nusselt number, as fol- 
lows (Bej an and Kraus[|2003| >: 

0.023/?e°- 8 Pr - 4 . 



Nu 



(10) 



Re and Pr are the Reynolds number and Prandtl num- 
ber of the fluid, respectively. We assume that these 
two numbers and the heat transfer coefficient are con- 
stant along the wellbore axis, and are determined by 
the fluid velocity and fluid properties evaluated at a 
temperature of 60 °C. 

The main interest of the present study is the stabil- 
ity of the wellbore. As explained in |Zoback| ( |2007| ), 
an unstable wellbore is defined as a borehole in 
which failure produces so much failed material from 
around the wellbore that it cannot be eliminated by 
mud circulation. To measure the instability, we cal- 
culate the cooling of the rock located one borehole 
diameter (i.e., 10 cm) into the rock reservoir (away 
from the borehole wall) and use this temperature de- 
crease to check the failure criterion in Eq. ([7]). If the 
failure criterion is fulfilled, we expect a failure zone 
with a maximum thickness of approximately 10 cm 
around the wellbore. In this case, the thermal frac- 
turing can potentially produce failed materials with 
sizes on the order of the borehole diameter. There- 
fore, as a borehole instability criterion, we check to 
see whether the failure criterion in Eq. ([7]) is fulfilled 
at a depth of one borehole diameter into the rock. 

To examine the convergence of the COMSOL pro- 
gram, we performed a simulation with a double- 
resolution grid and found that the difference between 
the two temperature solutions was less than 1%. 

3.2. In Situ Stress State 

Consider a reservoir under the large horizontal tec- 
tonic stresses typical in the studies of Fennoscan- 
dia ( |Stephansson et al.[|1989| ). The maximum (cr H ) 
and minimum (<r h ) principal horizontal stresses are 
given by ( [Fejerskov et al.[ |2000[ |Stephansson et al. 
T9891 ): 



cr H (z) = 2.8 + 0.04z MPa, 
cr h (z) = 2.2 + 0.024z MPa. 



(11) 
(12) 



The maximum horizontal stress is oriented north- 
west/southeast (see the red, dotted line in Fig. [3]). 
The vertical stress is given by <T h /<T V =1.1. 



We assume that the water in the wellbore diffuses 
into the low-permeability rock surrounding the bore- 
holes. Bernoulli's principle states that the velocity 
of the fluid results in a pressure decrease. The ve- 
locity yields a relative correction to the hydrostatic 
pressure (on the order of v 2 /gz). At the depths (i.e., 
greater than 3000 m) and fluid velocities (on the or- 
der 1 m/s) considered in this paper, the relative ve- 
locity correction to the fluid pressure is on the order 
10" 4 . We, therefore, neglect the difference between 
the borehole fluid pressure and pore pressure in the 
rock, i.e., Ap = 0, and set the pore pressure in the ef- 
fective stress tensor equal to the hydrostatic pressure, 
i.e., P p = p/gz, where g is the gravitational accelera- 
tion. 

We use values for the Young's modulus and 
Poisson's ratio that are typical for Fennoscan- 
dia ([Stepha nsson et alTj [1989) and a thermal ex- 
pansion coefficient that corresponds to a silica con- 

10 4 



tent of 73% (Zoback 

a t 



2007): E = 6.0 



MPa, 



8.5 • 10" 6 K" 1 , and v = 0.2. For simplicity, the 
Biot- Willis parameter is set to one (b = 1.0). 

4. Results 

First, consider the minimum cooling of the bore- 
hole wall required to induce thermal tensile fractures. 
When the wellbore is oriented along one of the prin- 
cipal stress directions and Ap f = AT = 0, the mini- 
mum value of the effective stress in Eq. ([7]) becomes 
o~ t min0mm) = 3<t 2 2 - o~w A wellbore that is directed 
along the cr h , cr H , and cr v directions, results in the 
minimum effective stresses: 



(T„ 



min 

V 



3<Xy - CTji 
3<T V - <T h 
3<T h - (Tjj 



(13) 
(14) 
(15) 



For the principal stress values and depths considered 
here, we find that cA n < cr^ in < cr^. Thus, we 
expect thermal failure to be most dominant in well- 
bores that are oriented along the minimum horizontal 
stress direction. This premise is verified by Fig. [3j 
which shows the minimum AT that induces thermal 
failures for all possible wellbore orientations. The 
results are calculated from Eq. ([7]) for three different 
depths: 3000 m, 4000 m, and 5000 m below the sur- 
face. The secant method is used to find the AT value 



6 



that yields cr t min = 0. We see that the magnitude of 
AT depends strongly on how the wellbore is oriented 
with respect to in situ stresses. As expected from the 
above arguments, the minimum AT appears for well- 
bores that are directed along the minimum principal 
stress direction. Wellbores that are tilted along the 
maximum stress direction require the largest values 
of AT to initiate thermal failure. The maximum AT, 
which is above 130 °C for depths greater than 3000 
m, is more than three times the minimum AT. The 
large AT values along the cr H direction imply that 
thermal failure is completely absent for wellbores 
oriented along this axis. In contrast, for wellbores 
oriented along the minimum stress direction, thermal 
fracturing is expected for temperature changes in the 
range of 40-60 °C. In deep geothermal systems with 
reservoir temperatures greater than 140 °C, a temper- 
ature change of this order is feasible. Thus, for the 
model system described in Section [3} we expect ther- 
mal fracturing to be present for boreholes tilted along 
the minimum stress direction. 

To study thermally induced wellbore instabilities, 
we performed a numerical COMSOL simulation of 
the reservoir's temperature evolution. Fig. [4] shows 
the minimum value of the wellbore instability crite- 
rion along the 2000m-long wellbore segment for all 
possible orientations of the wellbore after one year 
and 20 years of geothermal production. In each fig- 
ure, the blue lines confine two blue pockets that rep- 
resent the critical orientations for which thermally 
induced instabilities are present. Within less than one 
year of geothermal production, borehole instabilities 
occur for the wellbores that are tilted by an angle of 
6 > 45° along the minimum stress direction. After 
20 years of production, this critical region expands 
to include wellbores that are tilted by an angle of ap- 
proximately 6 > 35°. Thus, for the most part, the set 
of critical borehole orientations is mainly fixed after 
one year of production, and only minor changes in 
the set are observed after 20 years of production be- 
cause the cooling of rock close to the borehole begins 
as a rapid process and quickly equilibrates towards 
a quasi- stationary process in which the temperature 
changes slowly. 

We also studied of how the instability criterion de- 
pends on the injection rate and injection temperature 
of the fluid. Fig. [5^ shows the time evolution of the 



120 




200 



130 



Figure 3: (a) All possible wellbore orientations are represented 
as points on the unit disk. A well that intersects the lower hemi- 
sphere of the unit sphere at the point x, which is parameterized 
by the angles and 0, is mapped to the point x by a projection 
onto the xy plane, (b) The temperature change of the rock (at 
the borehole wall) required to induce thermal fractures for all 
possible orientations of the wellbore at a depth of 3000 m. (c) 
The temperature change of the rock (at the borehole wall) re- 
quired to induce thermal fractures for all possible orientations 
of the wellbore at a depth of 4000 m. (d) The temperature 
change of the rock (at the borehole wall) required to induce 
thermal fractures for all possible orientations of the wellbore at 
a depth of 5000 m. In (b) - (d), the red dotted line shows the di- 
rection of the maximal horizontal stress, while the inner (outer) 
black circle represents wellbores that are tilted by an angle of 
30° (60°) with respect to the z axis. 



minimum value of the wellbore instability criterion 
along the 2000m-long borehole segment for injection 
rates of 5, 15, and 30 L/s. The wellbore is oriented 
along the direction 6 = 60°, (j) = 75° and the in- 
jection temperature is 45 °C. The figure shows that 
the value of the stability criterion increases weakly 
with decreasing injection rate. This result is logi- 
cal, given that a lower injection rate implies a lower 
rate of heat extraction from the reservoir. The phys- 
ical reason for this weak injection-rate dependency 
is that the heat extraction rate is predominantly gov- 
erned by the heat transfer coefficient, h, which scales 
with the Reynolds number of the fluid as h oc Re *. 
Fig. [5J3 shows the time evolution of the minimum 
value of the wellbore instability criterion along the 
2000m-long borehole segment for injection temper- 
atures of 30 °C, 45 °C, and 60 °C. The orientation 




Figure 4: (a) The figure shows the minimum value (along the 
well segment) of the wellbore instability criterion for all possi- 
ble orientations of the wellbore after one year of operation, (b) 
The figure shows the minimum value (along the well segment) 
of the wellbore instability criterion for all possible orientations 
of the wellbore after 20 years of operation. In both plots, the red 
dotted line shows the maximum horizontal stress direction. In 
each plot, the two pockets that are confined by the blue dotted 
lines represent the orientations for which we expect thermally 
induced borehole instabilities. The injection temperature of the 
water is 45 °C, and the injection rate is 15 L/s. 



of the wellbore is the same as in Fig. [5^, and the in- 
jection rate is 15 L/s. We see that the value of the 
instability criterion depends strongly on the injec- 
tion temperature of the fluid. The fluid-temperature 
dependency is mainly determined by the elasticity 
parameters of the reservoir. For the model system 
studied in the present paper, a temperature change 
of 15 K yields a thermal stress change on the order 
of a t El5/(l - v) = 9.6 MPa. This estimate is in 
good agreement with the injection-temperature de- 
pendency shown in Fig.[5]3. As a consequence of this 
large thermal stress-change, the set of critical bore- 
hole orientations is strongly affected by the injection 
temperature of the fluid. Adding or subtracting a 
thermal stress of ~ 9 MPa in Fig. [4] provides an esti- 
mate for the case where the injection temperature is 
60 °C or 30 °C. Even when the injection temperature 
is as high as 60 °C, a significant number of borehole 
orientations become unstable. 

Our numerical simulations reveal two import con- 
siderations regarding long-term thermal effects on 
EGS wellbores. First, the system is extremely sen- 
sitive to how the wellbores are drilled with respect 
to the in situ stress state. For boreholes that are ori- 
ented along one of the critical directions, thermally 
induced fractures with risk of causing wellbore in- 
stabilities are observed after only a short time of pro- 
duction via the geothermal system. In other words, 
prior knowledge of the reservoir's stress state is cru- 
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Figure 5: (a) The figure shows the time dependence of the 
minimum value of the wellbore instability criterion along the 
wellbore for injection rates of 5, 15, and 30 L/s. The injection 
temperature of the water is 45 °C. (b) The figure shows the time 
dependence of the minimum value of the wellbore instability 
criterion along the wellbore for injection temperatures of 30, 
45, and 60 °C. The injection rate is 15 L/s. In both plots, the 
wellbore has an orientation of 6 = 60°, 6 = 75°. 



cial for optimal design. Second, the extent to which 
thermal fracturing appears is largely determined by 
the temperature of the injected fluid. This strong 
injection-temperature dependency opens a path for 
controlling the thermal failure process. Our results 
provide important insights regarding long-term ther- 
mal effects that are essential for optimal engineer- 
ing of WEGSs as well as geothermal wells and oil 
wells, in general. Importantly, our results can also 
be interpreted in terms of CCS technology. One 
of the major concerns regarding CCS technology is 
fracturing of the caprock caused by the C0 2 injec- 
tion ( [Zoback and Gorelick[ |2012| ). Recent studies 
have indicated that the thermal stresses induced by 
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injecting cold C0 2 play an important role in deter- 



mining whether the caprock will fracture ( |Luo and 
Bryantj [20T0t |Goodarzi etaTj |20T0l |20T2| ). Our re- 



sults imply that drilling through the caprock along an 
optimal direction reduces the risk of thermal caprock 
fracturing close to the well. 

5. Conclusions 

This paper contains a theoretical study of thermal 
fracturing and considers the effects of borehole ori- 
entation in deviated geothermal wellbores. We found 
that the risk for thermally induced wellbore instabil- 
ities is highly sensitive to how the wellbores are ori- 
ented with respect to the in situ stress state of the 
rock reservoir. Thermal failures can be completely 
eliminated if the system is optimized with respect 
to the principal stress directions. In contrast, if the 
system is not optimized, thermally induced fractures 
with risk of causing instabilities are observed after 
less than one year of geothermal production. In ad- 
dition, we found that the degree of thermal fracturing 
depends strongly on the temperature of the injected 
water, while it is only weakly dependent on the in- 
jection rate of fluid. This strong fluid-temperature 
dependency provides a direct way to control the de- 
gree of thermal fracturing for a given reservoir stress 
state. 
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